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ABSTRACT 

We have performed a linear pulsational stability survey of 6 series of 
long period variable models with M = 1.0M Q , L = 3000 — 8000L©, and 
(X,Z) = (0.700, 0.020), (0.735, 0.005). The dynamic and thermodynamic couplings 
between convection and oscillations are treated by using a statistical theory of nonlocal 
and time-dependent convection. The results show that the fundamental and all the 
low overtones are always pulsationally unstable for the low-temperature models when 
the coupling between convection and oscillations is ignored. When the coupling is 
considered, there is indeed a "Mira" pulsationally instability region outside of the 
Cepheid instability strip on the H-R diagram. The coolest models near the Hayashi 
track are pulsationally stable. Towards high temperature the fundamental mode first 
becomes unstable, and then the first overtone. Some one of the 2nd -4th overtone 
may become unstable for the hotter models. All the modes higher than 4th (n > 4) 
are pulsationally stable. The position and the width of such an instability region on 
the H-R diagram critically depends on the mass, luminosity and metal abundance of 
the star. The overall properties of the dependence are: 1) For the same mass and 
luminosity the instability region becomes slightly wider and moves to lower effective 
temperatures as the metal abundance increases; 2) For a given chemical abundance, 
the instability region becomes wider and moves to the lower effective temperature as 
their luminosity increases or their mass decreases 

For the luminous red variables seated outside the instability strip the dynamic 
coupling between convection and oscillations balances or may even overtake the 
thermodynamic coupling. Turbulent viscosity can no longer be ignored for the 
pulsational instability of the low-temperature red variables. The effect of turbulent 
viscosity becomes more and more important for higher modes, and may finally become 
the main damping mechanism of the pulsation. 

Subject headings: Convection - stars: Oscillations - stars: variables: long-period 
variables 
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1. Introduction 

Observationally, there exists a group of low-temperature luminous pulsating red variables to 
the right of the Cepheid instability strip on the H-R diagram. They are the most heterogeneous 
red giants and super-giants belonging to both population I and II. On GOVS (General Catalogue 
of Variable Stars), these luminous red variables are categorized into 3 different types according to 
their variability: Miras (M), semi-regulars (SR) and slow irregular variables (L), among which, 
the first type (Miras) have been well studied. It is currently believed that the Miras are stars on 
the AGB stage of evolution. Detailed summaries of the observational properties and theoretical 
work on Miras have been given by Whitelock (1990) and Wood (1990a,b). 

The nonlinear pulsation of Mira variables has been considered by Keeley (1970), Rose and 
Smith (1972), Wood (1974), Tuchman, Sack and Barkat (1978,1979), Hill and Willson (1979) 
Bowen (1988) and Perl and Tuchman (1990). The linear analysis is still useful for the stability 
survey. The linear pulsation of Mira variables were thoroughly studied by many authors (Kamijo 
1962, Langer 1971, Fox and Wood 1982, Ostlie and Cox 1986, Balmforth Gough and Merryfield 
1990, Cox and Ostlie 1991, Gong, Li and Huang 1995), among which Ostlie and Cox (1986) and 
Gong et al (1995) completely ignored the coupling between convection and oscillations. The 
coupling theory used by Kamijo, Langer, Fox and Wood, and Cox and Ostlie is over-simplified. 
Balmforth et al had used Gough's local time-dependent mixing length theory of convection in 
dealing with the coupling. 

The red pulsating variables possess very extended convective envelopes. The convective energy 
transport in the H and He ionization regions well exceeds 99% of the total value. Convection 
overwhelms the ft-mechanism of radiation and becomes the principal excitation (damping) 
mechanism of pulsation. The local time-dependent theory of convection has been used to interpret 
the red edge of the Cepheid instability strip (Baker and Gough 1979, Xiong 1980). When the 
surface temperature decreases, the dynamical coupling (through turbulent pressure and turbulent 
viscosity) between convection and oscillations becomes more and more important (Xiong 1977, 
Gough 1977, Stellingwerf 1984). For red stars outside the Cepheid instability strip, the dynamical 
coupling between convection and oscillations becomes as powerful as the thermodynamic coupling. 
Precisely speaking, the local time-dependent theory of convection can no longer (at least in a 
self-consistent way) be used to treat the dynamical coupling. For this reason, we have developed a 
nonlocal time-dependent statistical theory of convection (Xiong 1989). In this paper we assumed 
that the convection is quasi-isotropic. Such an assumption excludes turbulent viscosity, which is 
anisotropic. We have developed a more precise version of the nonlocal time-dependent statistical 
theory of convection, in which the dynamic equations of the third order correlation functions are 
derived, and the anisotropy of turbulent convection is also considered (Xiong, Cheng & Deng 
1997a). We have derived an equation of turbulent viscosity which is very similar to the Stokes 
equation of viscous fluid. 

Comparing with previous works, the present version stands on a solid base of hydrodynamics 
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for describing the dynamic behavior of convective motion. Therefore, it would provide a better 
approach to both the dynamic and thermodynamic couplings between convection and oscillations. 
Using this new theory, in this paper we have performed a linear pulsational stability analysis 
for the luminous red stars outside the Cepheid instability strip. The working equations of linear 
non-adiabatic oscillations are given in Sect Section |3| describes numerical results and theoretical 
arguments. A summary and concise conclusions of the present work are given at the end. 



2. Working Equations 

We have given a complete presentation of the nonlocal time-dependent theory of convection 
and all its mathematical formalism, i.e. the radiation-hydrodynamic equations of stellar 
oscillations (Xiong, Cheng & Deng 1997a). Among all the 14 equations, 3 deal with fluid motion 
(the conservation laws of mass, momentum and energy), 2 describe the radiation field (radiative 
transfer and energy conservation of the radiation field), and the other 9 equations are given for 
the second and the third correlation functions of turbulent velocity and temperature for turbulent 
convective motion. Within the frame work of the Eddington approximation, the departure from 
radiative equilibrium has been accurately taken into account. The gas and radiation are treated 
separately, and they become coupled through the absorption and emission processes of matter. 
In the mean time, careful consideration has been given to the anisotropy of turbulence and the 
inertia of turbulent convection. The turbulent viscosity, which is nearly identical in form to 
the Stokes formula for viscous fluid, is included automatically in the dynamic equation of the 
second-order correlation of turbulent velocity. When we adopt a simplified gradient-type diffusion 
approximation for the third-order correlations, this reduces the number of radiation-hydrodynamic 
equations down to 11. We have computed the linear non-adiabatic radial oscillations of p-modes 
for a solar model using these simplified radiation-hydrodynamic equations. The results show that 
the effects of departure from radiative equilibrium are negligible only except for very high order 
p-mode (n > 25) (Cheng & Xiong 1997). In this work, we will ignore the departure from radiative 
equilibrium. In the context of the present study this is safe, because the giants and super-giants 
almost always pulsate at low-order modes. Assuming that the gas and the radiation field are 
always in equilibrium, we can treat them together. In this way, the number of equations is reduced 
down to 10 and the pulsational equations are greatly simplified. Under such circumstances, the 
linear non-adiabatic equations can be written as 
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are the relative complex amplitudes of pulsation, x, Z and V are the auto- and cross-correlation 
functions of turbulent velocity and temperature fluctuation respectively, 
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and r* = jqT c . t c is the inertial time scale for turbulent convection, which is 
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while Li, L3 and L5 are the convective variables relevant to the third correlations. c\ and 
C2 are respectively the two convective parameters relevant to the dissipation and diffusion of 
turbulence. In our statistical theory of turbulent convection (Xiong 1980,1989), l e = C\Hp (Hp 
is the local pressure scale height) is the average size of energy-containing eddies, e = 2r] e x 3 /l e 
(r] e is the Heisenberg eddy-coupling constant) is the turbulent disspation (Hinze 1975). / = C2Hp 
is the turbulent diffusion length scale. The e-folding length of convective overshooting is about 
\A^/c\CiH-p (Xiong 1985). a and b are the coefficients of expansion and compression for the gas. 
Variables with a subscript "P" (or "T") denote the logarithm partial derivative with respect to 



More detailed definitions concerning the above expressions have been given in our previous work 
(Cheng & Xiong 1997). One thing that should be made clear is that the gas and radiation are 
treated as a whole (P = P g + P r ) in the present work, differing from our previous manner. 
Therefore, the contribution by radiation to all the material properties and to all the thermodynamic 
quantities, such as S, a, C p and V a d, has been included. 

The boundary conditions at the bottom layer are: 
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and the boundary conditions at the surface are, 



2/3 = 1, 



yi + U + 



r 3 w 2 



2/3 = 0, 



4y2 + 2y 3 -2/4 = 0, 



2/6 = \/ „ , ~ ^ (2 + iUT c )y 5 , 



1 + 2^/3% 



(18) 
(19) 
(20) 
(21) 



2/s 



12^3% 



(22) 



2/io 



12^3r? e 



2 + 



\ 1 + 2V3r? e V V3r?, 



+ wr c h/g, 



(23) 



As we have done in the previous paper (Cheng & Xiong 1997), the lower boundary is set 
in the convective overshooting zone, and the upper one is taken to be optically thin enough in 
the stellar atmosphere. Explicitly, we have adopted an optical depth of r = 0.01 for the upper 
boundary in the present work. 

The convective boundary conditions Eqs. (ffq)— (ff^) and (pl|)-(p3|) are given by the asymptotic 
analysis of the convective quantities within the overshooting zone. For their details we refer to 
our another work by Xiong et al (1997b), where the spatial oscillations of convective variables are 
amply discussed. 



3. The Numerical Results 

We have calculated the linear non-adiabatic pulsation for 6 series of models of luminous red 
stars. Their masses, luminosities, chemical compositions and of pulsational instablity region are 
listed in Table 1. The convective parameters are c\ = C2 = 0.75, this leads about the same energy 
transport efficiency as the original Vitense theory (Vitense 1958) when the mixing-length is taken 
to be one and half of the local pressure scale height, i.e. I = 1.5Hp, and the e-folding length of 
convective overshooting is about lA^/cic^Hp ~ 1.05Hp. The reddest model in the series is very 
close to the corresponding Hayashi track, while for the hottest model the convective flux is much 
smaller than the radiation flux in the second ionization region of helium. For these hotter stars, 
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the K-mechanism functioning in the second ionization region of helium has already operated, and it 
has become the main excitation mechanism. Therefore, our models actually cover all the possible 
temperatures of the luminous red stars outside the Cepheid instability strip, for the present setting 
of mass, luminosity and chemical composition. 

In the present work we use a simplified MHD equation of state (Hummer and Mihalas 
1988, Mihalas et al. 1988, Dappen et al. 1988). The neutral helium has been considered as a 
hydrogen-like atom for convenience of calculation of its energy levels. A significant departure only 
appears in the ground and the lower excited states. An analytic approach to the OPAL tabular 
opacities (Rogers and Iglesias 1992) and the low-temperature tabular opacities (Alexander 1975) 
is used for the calculation of opacity. 



3.1. "Mira" Instability Strip 

Table 2 gives the linear amplitude growth rates of the model series 2 of luminous red stars, 
r\ = —2TTLVi/cj r , where u> = uj r + iuoi is the complex angular frequency of linear non-adiabatic 
oscillations. The first column is the serial number of the model, the second is the effective 
temperature T e , the 3rd- 7th columns are the amplitude growth rates for the fundamental through 
4th overtones (taking into account the coupling between convection and oscillations). We have 
actually calculated the first 12 modes (n = — 11) for non-adiabatic oscillations. However, all 
the modes with n > 4 are pulsationally stable for luminous red stars. For the sake of clarity, 
we give in Table 2 only the first five modes. Aimed at examining the effect of convection on the 
pulsational instability for red stars, we have also computed the non-adiabatic pulsation ignoring 
the coupling between convection and oscillations (however, convection has been included in 
the static models). If convection does not vary during the oscillation of a star, the convective 
variables y§ = y% = y-j = ys = yg = yw = 0, Eqs. (|])-(|l0|) and the convective boundary conditions 
Eqs. (p^)-([T7|) and (|21~|)-([23|) can be omitted. Omitting the terms related to convection, Eqs. (||) 
and ([l^) are reduced to 

(P yi )-(4^P+prW)y 3 = 0, (24) 
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The number of pulsational equations reduces down to 4. That is the decoupling between 
convection and oscillations. Eqs. ([l]), (J3j) (|24|), fl25|) and the boundary conditions Eqs. (13), (14), 



and (|18[) (J20|) were used to calculate the non-adiabatic pulsation of stars ignoring the coupling 
between convection and oscillations. The corresponding growth rates of the first five modes 
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(n = — 4) are given as the 8-1 2th columns in Table 2. Depicted in Fig. |T| are the variations of 
the amplitude growth rates of the fundamental mode and the first overtone against the stellar 
effective temperatures for the 6 model series. 

It is clear from Table 2 and Fig. [I] that, when excluding the coupling between convection 
and oscillations, the fundamental mode, the first overtone and most of the low order overtone are 
pulsationally unstable for all the models of luminous red stars. But this is not true. Obviously, 
this is due to the ignorance of the coupling between convection and oscillations. More generally 
speaking, no good interpretation can be made for the red edge of the Cepheid instability strip if 
one does not take such coupling into consideration (Iben 1966, Xiong 1980). 

When the coupling between convection and oscillations is considered, most of the low 
order modes and all the higher order modes (n > 4) become pulsationally stable. However, 
there does exist a pulsationally unstable region for the luminous red stars outside the Cepheid 
instability strip. The coolest models near the Hayashi track are pulsationally stable. Towards 
high temperature the fundamental mode first becomes unstable, and then the first overtone. Some 
one of 2nd -4th overtones may become unstable for the hotter models. All the modes higher than 
the 4th (n > 4) are pulsationally stable. For simplicity, and to make it distinct with the other 
pulsationally unstable region, we will call it "Mira instability strip" for the moment, though such 
a name may not be very appropriate. This instability strip includes not only the regulars Miras, 
but also the semi-regular and irregular variables. 

Our calculations show that the pulsational instability of a luminous red star depends critically 
on stellar metal abundance, luminosity and mass. The overall properties are: 

1. For the same mass and luminosity, the instability strip becomes slightly wider and its 
location moves to lower effective temperature in the H-R diagram as the metal abundance 
increases. 

2. For the same abundance, the instability strip becomes slightly wider and its location moves 
to lower effective temperature in the H-R diagram as the the luminosity increases or the 
mass decreases. 



The above characteristics of luminous red variables are closely connected with the structure 
of the convective zones of these stars. Convection is the dominant factor that controls the 
pulsational instability for red stars. H2O and TiO contribute most of the absorption in the 
low temperature region, therefore the low temperature opacity strongly depends on the metal 
abundance. Moreover, the structure of the stellar convective zone is determined by the nature of 
the outermost super-adiabatic layer, hence the extension of the convective zone depends critically 
on the metal abundance. That explains why the Mira instability strip and the Hayashi track 
move to lower temperature as metal abundance increases. In the mean time, the structure of the 
convective zone depends also on the mass and luminosity. The above reasons make it easy to 
understand the properties of the Mira instability strip. 
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3.2. Excitation and Damping of Oscillations 

For the luminous red stars outside the instability strip, the convective energy transport 
exceeds 99% of the total quantity in the ionization regions of hydrogen and helium. Fig. ^ 
and Fig. [| give the integrated work Wp versus depth for the luminous stars pulsating at the 
fundamental mode. In Fig. ^, convection has been included for the static model, while the coupling 
between convection and oscillations is ignored in calculations of non-adiabatic oscillations. The 
only excitation and damping mechanism results from radiation. It can be found from Fig. ^ that 
for the two cooler models the excitation comes mainly from the the outermost gradient region 
of radiation flux, where x T approaches its maximum. This excitation mechanism functioning 
in the gradient region of radiation flux will be detailed in our second paper (Xiong, Cheng & 
Deng 1997b). The radiative excitation arising in the second ionization region of helium becomes 
negligible, because the radiative flux is far smaller than the convective flux here. The radiative 
excitation occurring in the second ionization region of helium becomes unnegligible for the hotter 
model. For a clearer display, the variation of \t an< ^ 1°§ {L r /L) m the stellar interior are given in 
Fig. |2] as well. Due to the convective energy transport, the radiative damping in the deep interior 
of a star is greatly weakened. Hence, the star is always pulsationally unstable when the coupling 
between convection and oscillations is ignired. 

Fig. ^ plots the quantities for the same luminous red star model, but with the coupling 
between convection and oscillations being considered. Wp, Wpt and W v i s represent respectively the 
contributions of the gas (and radiation) pressure, turbulent pressure, and turbulent viscosity. Their 
exact definitions are referred to our previous work (Cheng & Xiong 1997). W a ii = W p + Wpt + W v i s 
is the total integrated work, the value of which at the surface of a star should be the amplitude 
growth rate r\ = —2iruJi/u> r given by the non-adiabatic pulsation calculation. Our results show 
that they are very closely matched indeed. 

By comparing Figs. || and |j[ one can see that their characteristics are totally different. It can 
be found from Fig. ||] that with consideration of the coupling between convection and oscillations 
the excitation (and damping) region of pulsation, where W a ii changes obviously, extends deeper 
into stellar interior in comparition with Fig. 2. This is surely a result of the coupling between 
convection and oscillations. 

Convection leads to energy and momentum exchange inside a star, and hence affects the 
pulsational instability for variable stars. We call the effects of convective energy transport on 
variable stars as thermodynamic coupling, while the effect of momentum exchange on variable 
stars as dynamic coupling. 

The gas (and radiation) pressure term Wp in the total integrated work contains the effect 
of energy transport not only by radiation, but also by convection (i.e. the contribution by 
thermodynamic coupling between convection and oscillations). It is hard to unambiguously 
separate the effects of convective energy transport from radiative transport, because the two 
factors are closely coupled in the process of stellar oscillations. Following the conservation law 
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of energy (momentum), any change in convective flux (turbulent pressure) must lead to the 
corresponding change in the radiation flux (gas pressure) in order to balance the previous factor, 
and visa versa. By the same reason, it is not so easy to separate unambiguously the contribution 
of gas (and radiation) pressure Wp from that of the turbulent pressure Wpt- 

The combination of Wpt and W v i s is the dynamic coupling between convection and 
oscillations. The effect of turbulent viscosity is to convert the kinetic energy of ordered pulsation 
into random turbulent kinetic energy, and this process happens in the low wave number region of 
turbulent energy spectra . Through the nonlinear effects of fluid dynamics, the energy gained in 
the pulsational motion is cascaded into higher and higher wave number regions, and is eventually 
converted into thermal energy at the smallest turbulent eddies. Therefore the turbulent viscosity 
always works as a damping factor to pulsation, i.e. W V i S is always negative. Although W v i s is 
far less than Wp and Wpt in absolute values for the fundamental mode of luminous red stars, 
we cannot ignore its effect on the stability of pulsation. The contribution of turbulent viscosity 
increases very quickly towards higher modes for luminous red variables because viscous dissipation 
is proportional to the square of velocity gradient which grows as the oscillation mode increases. 
The turburlent viscosity becomes the main damping mechanism of the high overtones. This is the 
reason why the luminous red stars become pulsationally stable at high modes. Fig. ^demonstrates 
the integrated work diagram for the first, 3rd and 5th modes of a T e = 2550-fT model. Fig. [|a and 
Fig. H] clearly show the quick increase of W v i s towards higher modes. 

The coupling between convection and oscillations depends critically on the ratio of the 
time-scale of convective motion t c to that of pulsation r p = 2ir/uj r , i.e. iv r r c . When r c is much 
smaller than r p , i.e. oj r r c <C 1, the convection is in nearly quasi-steady variation. When lo t t c 
increases, convection lags behind the pulsational variation. The result is that not only its 
variation amplitude decreases, but also a phase lag exists. The interaction between convection 
and oscillations depends on both the amplitude of the convective variations and the phase of 
such variations. As we will show later on (Xiong, Cheng & Deng 1997b), the r c distribution 
in the interior of luminous red giants is very different from that of the RR Lyrae and the 
horizontal branch red stars. This difference results from the difference of their mass- luminosity 
ratios between them. The long-period variables have very low mass-luminosity ratio. Therefore, 
their internal structure has very high central concentration. This may explain the difference in 
pulsational characteristics of these two kinds of stars. In Figs. ||| and || we have also given the 
curves for uj r T c versus depth (logT). This shows that, for such luminous red stars, u> r T c < 1 holds 
for most of the interior of envelope model except the extremely outermost layer. 



4. Conclusions and Discussion 

By using a nonlocal and time-dependent theory of convection, we have very carefully treated 
the coupling between convection and oscillations. The linear stability analysis for luminous red 
stars gives us the following results: 
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1. The coupling between convection and oscillations is the dominant factor for the pulsational 
instability. When this coupling is not considered, all the low-temperature red stars are 
pulsationally unstable, while when it is taken into account, a Mira instability strip will show 
up outside the Cepheid instability strip; 

2. Except for some slightly hotter models which may have pulsationally unstable second to 
fourth overtone modes, luminous red variable stars pulsate in the fundamental mode or the 
first overtone. All the modes higher than 4 are pulsationally stable; 

3. For the low-temperature red stars, the dynamic coupling between convection and oscillations 
is in the same order of magnitude as the thermodynamic coupling, and may even overtake 
the later; 

4. The effect of turbulent viscosity grows very quickly towards high overtones for the luminous 
red stars. For high overtones, it becomes the main damping mechanism. 

The troublesome spatial oscillations of the thermodynamic quantities are triggered when the 
local time-dependent convection theory is used in calculating the stellar oscillations. (Keeley 
1977, Baker <fe Gough 1979, Gonczi & Osaki 1980). In the present calculations of the nonlocal 
convection time-dependent theory the spatial oscillations still exist but are effectively controlled. 
For the model of luminous red variable stars, the spatial oscillations happen at the extreme outer 
layer of the atmosphere, where ujt c S> 1. In the stellar interiors, where ojt c < 1, the troublesome 
spatial oscillations do not appear. It is difficult to trace the spatial oscillations in Figs. || and 
||, since the convective energy transport is far less than the radiative energy transport, and the 
turbulent pressure is far less than the gas pressure in the atmosphere of a luminous red giant. It 
seems true that the spatial oscillations do not affect the pulsation instability of luminous red stars 
to a considerable level. 

As for the long disputed pulsational mode for O Ceti, some authors think that O Ceti is 
pulsating in the first overtone (Kamijo 1962, Keeley 1970, Wood 1974 & 1981, Tuchman, Sack 
& Barkat 1978, Tuchman 1991), while some others believe that it oscillates in the fundamental 
mode (Hill & Willson 1979, Willson 1982, Bowen 1988). The theoretical pulsation period P and 
pulsation constant Q for a luminous red star model, with L = 5000^0, X = 0.700, Z = 0.020, are 
given m Table 3. Fig. | gives the P and Q of the fundamental and first modes versus the effective 
temperature for 6 model series. Our theoretical value of the pulsation constant Q is very close 
to that of Ostlie & Cox (1986), but our Qi for the first overtone is systematically smaller than 
that of Wood (1982) especially at the long period end. Following our numerical results, there is 
no model having its first overtone Q\ that exceeds 0.045 days. Normally, Qi < 0.041. Assuming 
for O Ceti, M « M , T e = 2900if (Wood 1990a,b) and that it pulsates at the first overtone, then 
we have L ~ 12500Lq or Mj> ; ~ —5.5. This would make it 0.7 m more luminous than the value of 
Mboi = —4.8 given by the P-L relation for Miras in the Large Magellanic Cloud (Glass et al 1987, 
Hughes & Wood 1990). This is unlikely to be real. Wood (1996b) proposed another picture for O 
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Ceti. If Mira's P-L relation depends on the metal abundance, then the discrepancy between the 
observed and the theoretical Q values could be removed. Our results of linear pulsational stability 
survey support Wood's guess. The instability strip moves indeed towards the low temperature 
as the metal abundance increases, as shown by Table 1 and Fig. 1. Our calculations show that 
if the metal abundance for LMC Miras are only one-half of that of O Ceti the red edge of the 
Mira instability strip in the Galaxy will be lowered by about 350K compared with the LMC's. 
Therefore the magnitude of O Ceti derived from the LMC P-L relation will be 0.6 m more luminous 
than it actually is. Let us suppose that the effective temperature of O Ceti is T e = 2900-ftT and 
Mbol = — 4.2, then the pulsation period of the fundamental mode given by the non-adiabatic 
calculation will be 330 days. Following our linear stability analysis, it appears more likely for O 
Ceti to pulsate in the fundamental mode than in the first overtone. 

This work is supported in part by the contract number 19573018 of National Natural Science 
Foundation of China (NSFC). 
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Firure Captions 

Fig. 1. — The amplitude growth rate of pulsation for the fundamental mode and the first overtone 
varies with the effective temperature, with conv. and without conv. denote the cases when the 
coupling between convection and oscillations are included and not included respectively. 

Fig. 2. — The variations of the integrated work W p (solid line), T2 (dotted), x T (dashed), 
and log (L r /L) with respect to depth (logT) for three models of model series 2 (M = 5OOOM0, 
L = 5000-Lq, X = 0.700, Z = 0.020). The coupling between convection and oscillations is ignored 
in the pulsation calculation, while convection is included for the corresponding static models. The 
three horizontal lines are the locations of the ionization regions of H and He. 

Fig. 3. — The integrated work vs. depth (logT) for the same models as in Fig. 2, however 
the coupling between convection and oscillations is included. Wp, Wpt and W v i s are respectively 
the components by gas (and radiation) pressure, turbulent pressure and turbulent viscosity. The 
variations of log (L r /L) and log (wr c ) are also given. 

Fig. 4. — The same as Fig. [3l but for first, third and fifth overtones, The effective temperature of 
corresponding static model is T e = 2550-fT. 

Fig. 5. — The variations of the pulsation periods P and the pulsation constants Q for the 
fundamental mode and the first overtone, against the effective temperature. The coupling between 
convection and oscillations is included. 
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